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ABSTRACT 



We report the first identification of the optical bands of the B 2 £ + -X 2 £ + system of AlO in the red supergiant VY CMa. In addition to 
TiO, VO, ScO, and YO, which were recognized in the optical spectrum of the star long time ago, AlO is another refractory molecule 
which displays strong emission bands in this peculiar star. Simulating the bands of AlO, we derive a rotational temperature of the 
circumstellar gas of r iot =700 K. By resolving individual rotational components of the bands, we derive the kinematical characteristics 
of the gas, finding that the emission is centred at the stellar radial velocity and its intrinsic width is 13.5 kms~' (full width at half 
maximum). It is the narrowest emission among all (thermal) features observed in VY CMa so far. The temperature and line widths 
suggest that the emission arises in gas located within ~20 stellar radii, where the outflow is still being accelerated. This result con- 
tradicts equilibrium-chemistry models which predict substantial AlO abundances only to within a few stellar radii. We argue that 
non-equilibrium models involving propagation of shocks are needed to explain the observations. 

Key words. Astrochemistry; Stars: winds, outflows; Circumstellar matter, Stars: individual: VY CMa; Line: identification 



1. Introduction 

VY Canis Majoris (VYCMa) is a peculiar red supergiant 
with a very high luminosity of 3-10 5 Lq, an effective tem- 
perature, T e ff, of 3200- 3500 K, and an enormo us size which 
may be up to 3000 R Q dWittkowski et alj|2012h . It is thought 
to be an evolved massi ve star with initial mass of 25+1OM0 
dWittkowski et al.ll2.012h. That ma ss was reduced to the current 
17+8 Mq ( Wit tkowski et alj|2012l) due to mass-loss phenomena 
that occurred at different stages of its evo lution. The curren t 
mass-loss rate of the order of 10 4 Mq yr _I dDanchi et al.ll 19941) 
is probably the highest in VYCMa's history and resulted in an 
extended (~7") r eflection and emission nebula that partiall y ob- 
scures the object dSmith et alJl200 it [Humphreys et"alll2005h . So 
far unknown remains the mechanism that elevates matter high 
enough above the stellar surface to cool it down to tempera- 
tures at which dust grains can condense. Once dust is formed, 
radiation pressure on grains can drive the powerful circumstellar 
outflow, but such a "standard" scenario alone does not explain 
all th e observed features of the VYCMa's outflow (Smith TeTaT] 
1200 lb . It has been suggested that the highly complex structure 
of the nebula may be related to h uge convective ce lls on the 
surface of the enorm ous star (e.g. ISmith et ail l200lt see also 
lSchwarzschildl[T97l . 

Due to its extreme parameters and the possibility that 
it will explode as a core collapse supernova, VYCMa has 
been a subject of considerable astrophysical interest. This 
has recently intensified due to detections of new molecules 
in its cold circumstellar environment. While the first ra- 
dio detections of o xides OH, H2O, S i O dat e back to late- 
1960s /early- 1970s dEliasson & Bartlettl Il969t iKnowles et ail 
1969; |Buhl et all 1 19741) . it is the recent developments in mil- 
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limeter and submillimeter astronomy that have allowed to dis- 
cover the rich molecular spectrum of VY CMa at radio wave- 
lengths. Line surveys (i.e. systematic observations of sub stan- 
tial frequency ranges) like thos e of IZiurvs et alJ (f2007) and 
iTenenbaum et alj yOlO) (see also Ro ver et al.ll2010h allowed to 
fully explore the radio spectrum yielding 23 different molecular 
species (plus isotopologues) presently known. To date, the radio 
spectrum of VY CMa is the finest example of molecular richness 
in oxygen-rich gas around an evolved massive star. 

Optical observers have realis ed the uniqueness of the spec- 
trum of VYCMa much sooner dJovlll942t lHerbiglll970l) than 
radio-astronomers and it is the unusual optical spectrum that 
have triggered much of the attention garnered by this object 
for more than a half a century. In addition to strong atomic 
emission lines, whose presence is a rare phenomenon by itself, 
VYCMa displays strong optical emission in molecular bands. 
Such molecular emission so far has been identified i n electronic 
bands of TiO, VO, ScO and YO dH viand et al.lll969UWallersteinl 
1197 lb lHerbiglll974t IWallersteinlfl986l) . It indicates the presence 
of warm (a few hundred K) circumstellar material close to the 
star. All the oxides that are seen in emission in VY CMa are also 
those refractory molecules which are most often observed in ab- 
sorption in the photospheres of cool stars. Another such refrac- 
tory molecule, AlO, has never been reported in emission (nor in 
absorption) in the optical spectrum of VY CMa. 

However, pure rotational emission from the AlO radical has 
been recently detected in V YCMa at millimeter wavelengths 
(Tenenbaum & Ziurvs 2009). Being a refractory and also a rel- 
atively abundant molecule, AlO is of special astrochemical in- 
terest. It proceeds the formation of AI2O3 which is thought to 
be the first molecule to condense in an outflo w of a late-type 
star like VY CMa (e.g. lSharp & Huebnerl 19901) . Its observations 
bear the potential of gaining our understanding of dust conden- 
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sation in outflows of cool stars, in particular of the formation of 
alumina dust. Due to an unresolved fine and hyperfine splitting 
of the pure rotation lines, the small number of detected tran- 
sition s, and the poor angular res olution of the available radio 
data (Tenenbaum & Ziurvsll2009l) , the origin of AlO emission 
in VY CMa remains poorly understood. A potential detection of 
AlO at shorter wavelengths (i.e. in the infrared or in the optical 
regimes) can help to better establish how the molecule is formed 
and what its role is in dust formation. 

In an archival spectrum of VYCMa from 2001, we found 
strong emission features which undoubtedly belong to the 
B 2 1, + -X 2 E + system of AlO. We describe the observations and 
the overall spectrum of VYCMa in Sects.[2]and[3] We charac- 
terise the AlO emission briefly in Sect.[3]and compare it to our 
simulations in Sect.|H where we also give basic characteristics 
of the emitting gas derived in our analysis. We discuss the dis- 
covery of AlO optical emission and its implications in Sect. [5] 

2. Observations and data reduction 

The spectrum of VYCMa was taken with the Ultraviolet 
and Visual Echelle Spectrograph (UVES) at the Very Large 
Telescope (Paranal Observatory, Chile) on 2001 April 12. Two 
exposures of duration 2 and 10 s were taken. The spectrum cov- 
ers three spectral ranges, i.e. 4010-5220A (blue), 6460-8315 A 
(red), and 8430-10300 A (far red). The slit width of 1" provided 
a spectral resolution of about R-45 000. We acquired the spec- 
tra through the ESO archive and reduced them using the UVES 
pipelinq^ and IRAF using standard procedures for echelle spec- 
troscopy. With the slit length of 10" and the slit position angle 
of 0°, the nebula of VYCMa filled a substantial part of the slit 
and a careful correction for the background was applied. The 
stellar spectrum was extracted within an aperture with a width 
of 2". With seeing of a/8 (full-width half maximum, FWHM) 
during the observations, this aperture encloses most of the flux 
of the central source (the cut off was at 5% of the peak) and 
avoids most of the extended nebula. In this study, we focus on 
the blue spectrum for which the wavelength calibration is reli- 
able to within ~0.4km s _1 (3 rms) or better. Because no spectro- 
photometric standard star was observed with the spectral setup, 
the spectrum of VY CMa could not be properly calibrated in 
flux. To roughly put the final spectrum on a flux s cale, we used 
the sp ectral energy distribution of V Y CMa from iMassev et alJ 
(2006), which is based on observations five years after those with 
UVES. As we checked at the visual light curve by American 
Association of Variable Star Observerfl within uncertainties of 
0.4 mag, VYCMa had the same brigh tness on the d ate of the 
UVES observations and those of Masse v et al.l ((2006). 

3. The optical spectrum of VYCMa 

The extraordinary appearance of the optical spectrum of 
VYCMa has been extensively described in t he literature (e.g. 
Hvla nd et alJll969t IWallersteinl fl97ll 119861: iHumnhrevs etall 
2005ir [n particular, a spectrum obtained a year before the UVES 
observations, at a similar resolution, and covering the range 
5175- 9140 A was described in detail by Wallerst ein & Gonzalez! 
d200lh . We find no significant differences in both spectra in the 
overlapping region, and therefore here outline only briefly the 
main spectral features. The photospheric (absorption) spectrum 

1 http://www.eso.org/sci/software/pipelines/ 

2 http://www.aavso.org/lcg 



is dominated by molecular bands of TiO as expected for an early 
to intermediate type M-type star. On top of this photospheric 
spectrum, VY CMa displays a great variety of emission features 
which include atomic lines and molecular bands. The atomic 
emission is dominated by neutral atoms with lines of Ba II be- 
ing the only known exception. Strongest is emission from res- 
onance lines of K I, Na I, Ca I, and Rb I that form classical P- 
Cyg profiles. Many of the weaker emission features (e.g. from 
Ti I, Cr I, V I, and Fe I), exhibit inverse P-Cyg profiles or appear 
as pure emission lines. The spectrum contains also broad (non- 
photospheric) atomic absorption lines, which may be partially 
filled with emission. 

The UVES spectrum contains strong and numerous emis- 
sion bands of TiO and at least one band of VO (at 7867 A) 
w hich all appear unchang e d whe n compared to the spectrum 
of IWallerstein & Gonzal ez (2001). Unfortunately, our spectrum 
does not cover the stronge st bands of ScO and YO, which were 
present in the spectrum of IWallerstein & Gonzalez! d2001l) . and 
there is no evidence for emission of these molecules within our 
coverage. 

The optical spectrum of VYCMa short-ward of about 
5100 A, which is partially covered by the UVES spectrum, has 
not been explored by observers at high resolution. Between 
around 4800 and 5600 A, we found strong emission features 
with a well recognizable 'wavy' pattern indicating a molecu- 
lar rotation spectrum. We have identified these features as be- 
longing to the B 2 I. + -X 2 I. + (blue-green) system of AlO. A 
careful inspection of the spectrum revealed the presence of the 
Au = — 1, 0, 1, 2 sequences. The other well known sequence 
Av = ~2 with the main head above 5335 A is not covered by 
the UV ES spectrum and is too weak t o be seen in the spec- 
trum of IWallerstein & Gonzalez! (l200ll) . The UVES spectrum 
also does not reach the A-X system of AlO, which is located 
above 1.05 fim and we do not know any published spectrum that 
covers this range. In Fig.[TJ we compare the sections of the spec- 
trum containing the AlO emission with our best simulation of the 
bands (see below). The Au=+2 sequence is hardly recognisable 
and essentially appears only in two band-heads, i.e. (u", o')=(2,0) 
and (3,1), that look similar to atomic emission features. Each of 
the lower overtones with Ai>=-1 , 0, 1 exhibits three to four unam- 
biguously recognizable heads (i.e. with u"<3 and i>'<4) and the 
molecular emission can be traced in individual rotational com- 
ponents even up to ~100 A from the main head. The head of the 
(0,0) band, which reaches a flux level of 2.5 times the local con- 
tinuum, is the most prominent emission feature in the blue part 
of the UVES spectrum. 

4. Analysis of the AlO B-X system in emission 

We simulated the emissi on in the B-X system of AlO using ana- 
logues procedures as in lKamihski et alJ d2010l) . We assumed that 
the emission comes from a volume of homogeneous gas with no 
velocity gradient and that upper levels are populated by absorp- 
tion of photons from a central source. The pumping to upper 
levels is less effective for optically thick lines, because of the re- 
duction of flux from the central source, and this effect was taken 
into account allowing a derivation of gas optical thickness in the 
radial direction. The emitting gas is parametrised by vibrational 
and rotational temperatures (their meanings are explained later), 
column density in the direction to the star, and micro-turbulence. 
A given value of column density determines the optical thick- 
ness of a given transition and, as described later, it is the optical- 
thickness effects that allow a determination of the column den- 
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Fig. 1. The simulation of the B-X system of AlO (red line) is compared to the spectrum of VY CMa (blue line). The simulation was 
obtained for r iot =700 K and 7\,;b=2200K. A rescaled spectrum of Betelgeuse is also shown (grey line) to illustrate the contribution 
of photospheric (and circumstellar) features that have not been included in the simulation. From top to bottom, the panels contain 
sections of the spectrum dominated by the Au=2, 1, 0, and -1 sequences. The positions of AlO band-heads are indicated with red 
markers and labels (giving v" and v'), while some of the strongest features in the stellar spectra are indicated by grey markings. 



sity. In our approach, best fit parameters do not depend directly 
on the flux calibration nor the dust extinction in the nebula. In 
order to compare the simulation results to the observations, the 
simulated spectrum was smoothed with a Gaussian, what may 
be interpreted as a convolution with a macro-turbulence profile. 

The molecular data were taken from a compilation of 
sources: the list of lines w as built from spectroscopic con- 
stants in Sakse na~et al.l (120081): the electro nic transition moment 
was adopted from Zenouda et al.l (1 19991): Fran ck-Condon fac- 
tors were taken from ICoxon & Naxakisl d 1985b ; rotatio nal line- 
streng th factors were computed numerically following Hougen 
( 1970) and verified by comparison to the analytical formula in 
Kovacs (1969). Because Al has only one stable isotope, all the 
analysis was performed for the single isotopologue 27 A1 16 0. 

Since we deal with a 2 £ + - 2 £ + transition (Hund's case b), 
each rotational level, designated by the quantum number N, is 
split into two sublevels, F\ and F%, with total angular momenta 
given by J=N+l/2 and J-N-l/2, respectively. This is known 
as spin- or p-doubling. The band is formed mainly by the dou- 
blet branches P\, P2 and Ri,Ri- There are also satellite branches 
P Q\2 and R Q2i , which were included in the simulations, but their 
intensity drops very rapidly with increasing N and they have 
no practical meaning for our analysis. The p-doubling in the P 



and R branches is taken into account in our simulations follow- 
ing Saksena et al . (2008). Although there is some controversy in 
the literature on the s plitting constants y for v" > (see e.g. 
Launila & Bergll201 ll for the most recent review), it does not 
affect our analysis. 

FigureQ] presents the results of our best simulation. The 
shapes of the heads, the heads-intensity ratios, and the shape 
of the overall rotational profile are very well reproduced by the 
simulation. Some minor deviations between observed and simu- 
lated spectrum are due to an imperfect spectrum normalization 
and a presence of photospheric and circumstellar features, which 
were not taken into account in the simulation. To illustrate the 
shape and positions of the contaminating features, we compare 
the spectrum of V Y CMa to a spectrum of another red super- 
giant (Betelgeuse=o'Ori, Ml la) in Fig.Q] This r eference spec- 
trum was extracted form the Paranal UVES Atlas (Bagnulo et al. 
2003). The match between the photospheric and circumstellar 
absorption features in both stars is very good except that some 
features are slightly stronger in Betelgeuse (e.g. HB) or filled 
with emission component^ For VY CMA, the analysed spec- 

In choosing the reference photospheric spectrum, we did not at- 
tempt to find a good match to the spectral type of V Y CMa. The general 
appearance of its photospheric spectrum is much better reproduced by 
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trum can be satisfactorily explained by the combination of a cold 
photosphere and the emission spectrum of AlO. We also note 
that from the comparison of observations to the simulation, it is 
apparent that the AlO emission is not affected by any underlying 
photospheric absorption of AlO. If VY CMa follows the trend 
observed in gianto which display AlO abs orption only at spec - 
tral types later than M4 (very weak at M5) (Kee nan et al.l fl969>. 



AlO B-X (0, 0) 

P(23) f(24) P(25) 



AlO absorption is negligibl e in VYCMa, w hich is of spec- 
tral type M4-M5 ( e.g. IWittkowski et alj|2012h or earlier (~M2, 
Mass ev et al.l2 006). Because the Au=0 sequence is strongest and 
least contaminated by the photospheric and circumstellar fea- 
tures, it was the primary subject of the further analysis. 

By obtaining a grid of simulations with variable input param- 
eters and comparing it to the observations, we were able to con- 
strain some physical quantities of the AlO emitting gas. The ro- 
tational temperature (r rot ), which corresponds to the population 
of rotational levels within the vibrational states, can be derived 
from the overall shape of the rotational contours. Particularly 
useful for putting constraints on T mt is the observed position 
of the local maximum of the emission in the (0,0) band (out- 
side the main bandhead), which occurs at 4654+2 A correspond- 
ing to A=21±2 (Fig.Q}. It implies T ml between 600 and 900 K. 
Additionally, for rotational temperatures higher than ~800 K, the 
emission components of the R branch are too strong at the band 
origin, and below ~500 K the emission of the P branches form- 
ing the tail of the band is too weak. Therefore the rotation tem- 
perature implied by the simulation of the (0,0) band is in the 
range 7' lot =700±200K. However, the shape and intensity of the 
(1,0) and (1,2) bands is better reproduced by simulations with a 
higher rotational temperature between 900 and HOOK. 

Optical-depth effects also have an influence on the shape 
of the rotational spectrum. Because they are particularly evi- 
dent in the intensity of the band-head with respect to the rest of 
the band, it is possible to decouple the optical thickness effects 
from those of varying r iot . Our best fit resulted in optical depths 
of T hea d=3.0, 1.0, and 0.14 in the heads of the (0,0), (1,1), and 
(2, 2) bands, respectively. These values are equivalent to setting 
the AlO total column density to jVaio=6-10 16 crrT 2 . Because our 
simulations do not contain any realistic velocity gradient, these 
values are upper limits on the actual optical depths and the col- 
umn density if the emission arises in an outflow. The shape of the 
rotational contour is also somewhat influenced by the assumed 
value of micro- and macro-turbulence, but we kept those param- 
eters constant at FWHM=1 kms -1 (equal to thermal broadening 
at -1000 K) and FWHM=6kms _1 , respectively. 

The vibrational temperature (r v ;b) is indicated by the rela- 
tive strengths of heads within each sequence. We found satis- 
factory results for a broad temperature range between 1600K 
and 2800 K. Particularly good fits to data were obtained when 
different vibrational temperatures (but still within the specified 
range) were assumed for different bands. The best simulation 
with 7' 1 - ot =700K and 7\,;b=2200K is shown in Figs.Q]and[3] 

By comparing the simulation to observations, we were 
able to derive the bulk velocity of the AlO gas. Using the 
cross-correlation method implemented in IRAF's task fxcor, 
we derived heliocentric radial velocity of Vh e ii =41 + 1 kms 
(or equivalently Vlsr=21±1 kms -1 with respect of the Local 
Standard of Rest). This velocity is confirmed when peaks of in- 



spectral types close to M4, e.g. that of the giant HD 11695 from the 
same spectral atlas. 

4 Showing AlO emission, VY CMa is, however, more similar to Mira 
variables, in which Al O intensity does no t correlate with spectral type 
and is highly variable (Kee nan et alJl969T) . 
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Fig. 2. Rotational components of the B-X (0,0) band of AlO as 
observed in V Y CMa (blue histogram) and in our simulation (red 
line). The individual lines are unresolved doublets of members 
of the Pi and P2 branches. Their peaks are designated by ro- 
tational quantum numbers. Much fainter rotational components 
of the R branch can be seen in the simulation (only two were 
marked in the figure). 
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Fig. 3. The P(23) rotational component of the B-X (0,0) band 
of AlO in the simulation (red line) and in the observed spec- 
trum of VYCMa (blue histogram). The velocity is given with 
respect to the average wavelength of the doublet, whose split- 
ting is 4.2kms~'. The positions of the doublet components are 
shown by the two dashed lines in the upper part of the figure. 
Lines in the lo wer part show the hyp erfine splitting in both com- 
ponents (from Launila & Berg 20lTJ). 



dividual rotational components are compared to their rest wave- 
lengths. Moreover, we do not see any systematic shift in peak po- 
sitions with increasing excitation energy of the rotational com- 
ponents. 

Because of the good separation of the rotational components 
and their low optical thickness at high N (usually t <k1), we 
can examine the shape of individual lines for the kinematical 
profile of the AlO-bearing gas. The measured FWHM of lines 
between P(20) and P(30) of the (0,0) band is 14.4+1.2kms-'. 



4 



Kaminski, Schmidt, & Menten: Optical AlO emission in VY CMa 



This value is however affected by several splitting effects. The 
nuclear hyperfine splitting due to nuclear spin of 7=5/2, which 
was not included in the sim ulation, is known to b e very small 
(about 1.5 km s" 1 at P(23): lLaunila & Bergll201 lh . i.e. smaller 
than our spectral resolution. Additionally, each rotational line is 
a blend of two components belonging the P\ and P2 branches 
separated due to the p-splitting, which is very well known for 
the (0,0) band. Its magnitude changes linearly with rotational 
number and is negligible at low N (e.g. 0.3kms~' at N-2), but 
significant at high rotational numbers (e.g. 6kms~' at N-33). 
Correcting the measured widths for p-doubling gives a typical 
FWHM of 13.5±1.5kms _1 . To illustrate the different splitting 
effects on the line profile we show them graphically in Fig.[3]for 
the P(23) line for whi ch the hyperfine structure is known from 
lLaunila & Berg j2011l) . 

A comparison of the individual profiles to the simulation in 
Fig. |2] reveals a small asymmetry of the observed profiles, i.e. 
the red wing is more extended (by ~2kms~'). Although there 
is some uncertainty in the molecular data which we used in the 
simulation, we believe this effect is real. Additionally, on the 
blue side of each rotational component, a very weak absorption 
can be seen (it is more apparent when the spectrum is compared 
to the simulation - see Figs.|2]-[3]l. This and the profile asymme- 
try suggest that each rotational component may have a P-Cyg 
type profile. Some absorption can be also seen next to each of 
the strong band-heads (Fig.Q]), what supports this interpretation. 
The absorption is however very weak and the P-Cyg shape is 
claimed here only as a tentative result. 

5. Discussion 

5.1. Interpretation of the emission properties 

The central velocity of the optical AlO emission 
(Vlsr— 21±1 kms -1 ) is consistent with the velocity of the 
AlO millimetre lines (25+6 km s -1 ) and is in excellent agree- 
ment with the stellar systemic velocity known from optical and 
radio/(sub-)millimetre observations (VLSR=21.0+1.5kms~ 1 ) 
(e.g. [Smith 2004, and references therein). Although, the central 
velocity of the AlO emission is much more accurately measured 
in the optical (where splitting effects are smaller and a large 
number of rotational transitions at high signal-to-noise ratio 
is observed simultaneously), the centres of optical emission 
lines may be red-shifted due to scattering on dust grains in 
the thick and expanding nebula (Ivan B lerkom & van Blerkoml 
119781: Romanik & Leung] Il98ll) : this effect might have been 
observed in some absorptio n and emission feature s from the 
outer nebula of VYCMa ([Humphrey s et al.l 120051) . However, 
since the position of radio and optical lines of AlO are consis- 
tent within the uncertainties, this scattering effect must be very 
small in our case (<3kms -1 ), i.e. certainly much smaller than 
suggested theoretically by Ivan Blerkom & v an Blerkoml (Il978l) 
for photospheric absorption lines in VY CMa. 

The outflow of VYCMa observed in millimetre emission 
lines has been ascribed to at least three distinct kinematic al com- 
ponents dTenen baum et al. 2010) with the fastest material mov- 
ing at expansion velocity of at least 40kms~ 1 producing broad 
lines with multiple peaks. In addition to those wide features, nar- 
row emission lines centered approximately at the stellar systemic 
velocity have been detected, with the narrowest emis sion known 
arisin g from NaCl with FWHM=16+3kms~ 1 dMilam et all 
12007b : TiO is another molecule w hich displays such n arrow 
features with FWHM=18+2kms - 1 dKamihski et al.ll2012l) . The 
narrowest of the pure atomic emission lines present in UVES 



spectrum have a FWHM of 22 km s . All these narrow lines are 
signatures of gas located within the inner outflow, perhaps even 
below the dust formation zone. The width of the AlO emission of 
13.5 km s~' is smallest among all detected (thermal) features in 
V Y CMa's circumstellar environment indicating that AlO probes 
deepest observed layers of the material leaving the photosphere. 
With the typical full-width at zero intensity of 27+2 km s~' cor- 
rected for the p-doubling, the AlO emission implies gas expan- 
sion velocity of V ex p=13.5+2.0kms~ 1 . 

It should be mentioned that maser lines of SiO (and to a 
lesser degree of H2O) are known to exhibit even narrower emis- 
sion components with widths down to lkms -1 . These maser 
lines undoubtedly arise in the inner, accelerating envelope of 
VYCMa (within tens of mas from the star), but they usually 
do not trace the whole outflow, as most of thermal lines do, but 
rather individual clouds within the outflow. Moreover, due to a 
complex excitation mechanism, the SiO masers are observed in 
a ring-like feature surrounding the photosphere of VY CMa (e.g. 
Zhang et al. 2012) showing complex dynamics. Hence, the half- 
widths of individual lines do not correspond directly to the ex- 
pansion velocity. 

The slight asymmetry observed in AlO line profiles, with 
the red wing being more extended, is reminis cent of what is 
obser ved in the (sub-)millimetre lines of TiO (Kaminski et al.l 
2012). This asymmetry may suggest that even the inner wind re- 
gion exhibits some degree of asymmetry. However, exactly the 
same type of asymmetry is expected to arise due to scattering in 
an expanding dusty nebu la (Ivan B lerkom & van Blerkom 1978; 
Rom anik & Leung1ll98ll) . In addition to the optical thickness 
and geometry of the scattering medium, this effect is strongly 
dependent on the scattering properties of dust and for some dust 
types produces profile asymmetry wit hout any apparent shift in 
line centres dRomanik & Leungl 19811) . which could be the case 
for the AlO lines. The origin of the asymmetry is not clear 
and requires modelling including a detailed treatment of scat- 
tering on dust. This is beyond the scope of this work. We note, 
however, that if the observed profiles are re-shaped by scatter- 
ing, they are also slightly broadened and the intrinsic width of 
the AlO emission is even smaller than quoted above. The pure 
rotational lines of AlO obse rved at millimetre wavelengths by 
iTenenbaum & Ziurvsl ((2009), for which scattering can be com- 
pletely neglected, are severely affected by the hyperfine splitting 
and modest signal-to-noise ratios which complicate the profile 
analysis and the asymmetry cannot be verified. However, the au- 
thors were able to constrain the line widths to 10-15 km s _1 what 
is consistent with our analysis of the optical profiles. 

The vibrational temperature we derived describes the pop- 
ulation of vibrational levels in the upper electronic state in the 
sense of the Boltzmann distribution. Because the excitation is 
most likely radiative (the electronic state B 2 I. is 29691 K above 
the ground level), this approach is a crude simplification. The 
population of vibrational levels in B 2 E depends on the popula- 
tion in the ground state, the spectrum of the external radiation 
field, and the strengths of absorption bands populating the vibra- 
tion levels. By taking Franck-Condon factors as band intensities 
and assuming a flat stellar spectrum, we calculated that a tem- 
perature between 1200 and 1800K is more representative as the 
vibrational temperature. Taking the large uncertainties and the 
devious meaning of TVib, it is of no practical value for the further 
discussion. 

On the other hand, the rotational temperature of 700 K is 
expected to be close to the gas kinetic temperature, T^ n , what 
makes it a valuable parameter. T he rotational temperature we 
derive is much above that found bv lTenenbaum & Ziurvsl d2009) 
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(T rot ~250 K) on the basis of a rotational-diagram analysis. Their 
value is however very uncertain as only three lines were mea- 
sured at low signal-to-noise ratio. Our value of T Iot is consistent 
with temperatures derived from analyses of other optical bands 
of refractory molecules observed in VYCMa in emission — 
Phill ips & Davis! d!987l) found 7' rot =600± 150K from an analy- 
sis of the TiO j' system in emission, while iHerbigl (I 1 9741) found 
7/ rot =370 and 790 K for the orange (A-X) band of ScO observed 
at two different dates. 

5.2. The origin of the AIO emission 

In the 2-dimensional UVES spectrum, the AIO emission is un- 
resolved meaning that the emission region is smaller than 078 
FWHM. The width of the AIO profile allows to constrain even 
better the size of the emitting region. By combining tangential 
(proper ) and radial motions observed in H2O masers at three 
epochs, Richa rds et alj 0998) found expansion of the envelope 
within 75 to 440 mas from the star. By applying their expansion 
model (updated to the stellar distance of 1.2kpc) we get a radius 
r=18+6R* (R*=2000R Q ) for the AlO-emitting region. 

The type of observed line profiles also provides strong con- 
straints on the dimension of the AIO emission region. If the gas 
was much closer to the star than the 18 R* found above, say a 
few R+, the optical rotational profiles should display pronounced 
absorption components arising in gas seen in front of the stellar 
disk. Combined with a velocity gradient (which is observed in 
this region in SiO maser lines), such configuration should pro- 
duce P-Cyg profiles with a high absorption-to-emission intensity 
ratio, i.e. much more pronounced than those we observe (if any 
abso rption is present a t all). 

iDecin et all (120061) modelled the envelope of VY CMa with 
a variable mass loss and derived thermal model of the enve- 
lope with a detailed treatment of heating and cooling mecha- 
nisms. Assuming that 7\i n ~ T I0t , we find that the tem perature of 
700+2 00 K corresponds to a distance of 20+5 R* in lDecin et alj 
(2006), which is consistent with the constraints found above. In 
the same model, gas of such temperature is located in the accel- 
erating outflow. Their velocity law applied to the AIO expansion 
velocity (y exp =13.5kms _1 ) indicates location within ~10R*. 
The hydrogen density at 10-25 R* in that model (correspond- 
ing to the mass-loss rate of 10~ 4 MQyr~') is within «h=10 7 - 
10 9 cirT 3 . Because critical densities for the observed rotational 
transitions of AIO are between 5 TO 5 cm -3 (for N « 1) and 
5T0 9 cirT 3 (for N * 20fl the population of rotational levels 
has a chance to be thermalized by collisions in this region (cf. 
IWoitke et al.ll9 99). justifying our assumptio n that 7kin~ r rot . We 
must however note that the model of De cin et a 1. (2006) should 
be treated with caution as its parameters were constrained by ob- 
servations of CO lines with the highest transition being 7=7-6 
(155 K above the ground, « cl ; t =4-10 5 cirT 1 ) which does not probe 
well the accelerating outflow. 

The above discussion suggests that the AIO emission comes 
from a region of a diameter of about 40 R* (or 074). This al- 
lows us to calculate the AIO mean abundance. From our simula- 
tions, we get the AIO total column density of Na\o=6- 10 16 cm" 2 . 
This value has a factor of a few uncertainty resulting from the 
uncertainties in the other simulation parameters. For the mass- 
loss rate of 10~ 4 Mq yr -1 , the envelope extending from 1 R* to 
20 R*, the mean molecular weight of 1 .3, and the mean velocity 



5 The critical densities for the AIO rotational transitions in the X 2 ?. 
state were calculated for the collisional cross-section cr= 10~ 15 cm 2 and 
gas thermal broadening of 1 km s -1 . 



of 6 km s -1 , the total hydrogen column is A r H=3T0 24 cm" 2 giving 
AIO mean abundance of [f]-\og{NAio /Nh) »-7.7. This value is 
uncertain by at least 1 dex. 

Chemical equilibrium models of circumstellar gas around 
evolved oxygen-rich stars predict high abundances of AIO only 
at relatively high temperatures. For instance, in th e chem ical- 
equilibrium model calculated bv iGail & Sedlmavrl ([1998), the 
AIO abundance is expected to be higher than [f]-—9 only at 
temperatures between 1100 and 2000 K. Below these temper- 
atures, the abundance of AIO drops steeply with the distance 
from the star (being converted to AI2O, AI2O2, and AI2O3) and 
at 800 K is already 6 dex below its peak value ([/ max ] ~-8.5). 
Including condensation of alumina species into such models re- 
sults in even stronger depletion of A IO at temperatures lower 
than 1800K (Sharp & HuebneifT9 90). According to these mod- 
els, therefore, the AIO gas should be located in a relatively thin 
layer a few R+ above the photosphere. This is strongly contra- 
dicted by our constraints on the size o f the AIO emission region 
(see also [Tenenbaum & Ziurvs 2009) suggesting that the heav- 
ier alumina compounds, e.g. AI2O3, are formed less effectively 
than the models predict. Because AI2O3 is expected to condense 
to corundum, which is often thought to be the first c ircumstel- 
lar condensate (Shar p & Huebnerfl990tlLoddersl2003h . this dis- 
agreement may signalize that we should revise the picture of how 
effectively (if at all) alumina dust is formed in oxygen-rich stars 
like VYCMa. 

That the chemical-equilibrium models fail to explain the alu- 
minium c hemistry in the envelope of V Y CMa has been already 
noted in iTenenbaum & Ziurvs! d200 9). They proposed that the 
discrepancy may be caused by "freezing out" of the equilib- 
rium abundances at a few R* or that shocks disrupt the pro- 
cess of molecule fo rmation at larger distances. "Freezing out" 
dMcCabe et aljfl979b occurs when the time scale of chemical 
reactions (T c h e m 00 (Kri)~ l , where K is the reaction rate) be- 
comes much longer than the dynamical scale of the outflow 
(i"dyn ~ r/V exp ). As long as chemical reactions are faster than 
the dynamical time scale, the molecular abundances are thought 
to be well described by equilibrium values, but - as the den- 
sity drops - at some point in the outflow the reactions be- 
come inefficient and abundances do not change any more be- 
ing fixed at a cert ain level corre sponding to some equilibrium 
values. Following iMcCabe et alj d 19791) . we estimate that such 
freezing-out should occur in the outflow of VY CMa at densi- 
ties of abo ut «H=10 10 cm~ 3 , which are expected at a distance 
of few R* dDecin et al. 2006). Such a "freezing out" scenario 
in VY CMa would be supported for AIO by its abundance being 
close to the maximal value predicted by equilibrium models (and 
indeed those maximal abundances correspond to distances of a 
few R*). Whether such process can really take place in VY CMa 
should be verified by constructing a detailed model which would 
take into account the variable conditions in the outflow. 

We favour the alternative scenario in which the observed 
over-abundance of AIO at larger radii is a result of shocks related 
to the stellar activity. There are strong suggestions that shocks 
are present in the close vicinity of V Y CMa: 

- It has been proposed that the complex nebula surround- 
ing VY CMa was produced by a violent and episodic pro- 
cess that most likely involves creat ion of shocks ( e.g. 
episodic mass-ejections advoc ated in lSmith et al.ll200Tl and 
iHumphrevs et alJ 12005. 2007). This proposal, however, is 
based on observations of the extended nebula, formed by at 
least lOOOyr of mass-loss history, and it is not sure whether 
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the same kind of activity is currently taking place near the 
stellar surface. 

- The presence of extend ed and variable emissio n and broad 
absorption of Ha (e.g. Humphre ys et al.l 12005) within the 
otherwise neutral nebula indicates the presence of a process 
energetic enough to ionise hydrogen (although some of the 
emission, especially in the outermost nebula, may arise in 
the nearby H II region). 

- The star is known as an irregular variable, with changes of 
the optical light reaching 2 mag. This variabi lity may be par- 
tially explained by (non-radial) pulsations dHumphrevs et alj 
2005), which would be a potential driving source of shocks 
near the stellar surface. 

As an extreme supergiant very close to its Hayashi 



limit dWittkowski et al.l |20 1 2|). V Y CMa has huge convec- 
tive cells dSchwarzschildl whic h are expected to 

severely affect its circumstell ar medium (ISmith et alJl200U 
iHumphrevs eTaT]l2005l l2007l) . 

Circumstellar AlO bands are rarely observed. Their pres- 
ence and erratic behaviour h as been long known in Mira 
variables (Keen an et alj fl969l) and was reported in an enig- 
matic object UEqu (Barnbaum et al .1 119961) . They have been 
also observed in some members of the recently recognized 
family of nova-like variables called 'red novae'. These in- 
clude V4332Sgr, whic h has the optical a n d NIR bands 
of A lO in emission dBaneriee et al.l 120031: iTvlenda et al.l 
120051) . and V838 Mon, which shows circumstellar absorption 



in those bands ( Kaminski et "all 120091: iTvienda et alJ l201ll 
IE vans et al.l2003HLvnch et al.l2004l) . The AlO emitters seem 
to share one common characteristics - their circumstellar en- 
vironments have been shocked by violent motions, e.g. large- 
amplitude pulsations in case of Miras, and explosions in case 
of red novae and U Equ. The presence of shocks would jus- 
tify VY CMa belonging to the group of AlO-emitters. 



It would be indeed very surprising if the star were not to de velop 
shocks near its surface. As shown by ICherchneffl (|2006) (for 
asymptotic giant branch stars), propagation of shocks through 
material above the stellar photosphere changes its chemical con- 
tent considerably with respect to what is expected in chemical 
equilibrium. Future models attempting to explain the chemistry 
in VY CMa's envelope should explore this non-equilibrium ap- 
proach. Since shocks are expected to produce transient phenom- 
ena, there is also a possibility of observational verification. Most 
of the previously described molecular bands seen in VY CM a 
in emission are variable (e.g. I Wallersteinl fT986i: lHerbidfl974l) . 
although it is not clear whether the absolute intensity or the in- 
tensity relative to the variable continuum changes. It is also not 
known what are the shortest time scales on which the variabil- 
ity occurs except that it may be shorter than years. To explore 
the possibility that these optical emission bands arise due to the 
presence of shocks, it is highly desirable to verify whether the 
AlO bands are variable too and whether there is any correlation 
with the variability of the other molecular bands and hydrogen 
recombination lines. 
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